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1. INTRODUCTION 

One of the drawbacks of the graving dependence on air travel is the 
increased iirpact of aircraft noise. Assessment and reduction of this impact 
requires identification of aircraft generated noise levels. To this end, 
the Noise Technology Branch of NASA/Aircraft Noise Reduction Division is 
developing and refining an aircraft noise prediction computer program. 

Noise estimination includes consideration not onlv of noise sources on the 
aircraft, but also oi the propagation path between source and receiver. 

One of the factors affecting the noise transmission path is shielding of 
one jet by another. The shielding jet, because of the high temperature and 
flow speed with respect to the umediate surroundings, acts as a partial 
barrier between the shielded iet and the receiver. The resultant alter- 
ation of the propagation path not only affects the overall aircraft noise 
level, but also indicates the possibility of iet engine installation as a 
means of aircraft noise control. 

It is the purpose of this paper to discuss the development of the 
analytical model to estimate the shielding of one jet by an adjacent iet in 
a twin jet configuration. 

The problem of reflection and transmission of sound by a moving medium 
has been addressed assuming a plane wave incident on a plane interface (1, 
2,3,4). Ray tracing techniques have been applied to two-dimensional jets 
C Sj and cold jets (7). 

In the current studv, the three-dimensional case is considered. The 
noise source is a discrete frequency point source at rest with respect to 




the jet axis. The shielding jet is assuned to be a cylinder of heated flow 
in which the temperature and flow velocity profiles are constant across the 
jet. 

The three-dimensional character of the noise source is required in 
order to investigate the effect on shielding of the orientation of the 
emitting jet with respect to the shielding let. The finite cross-section 
of the shielding jet permits investigation of diffraction around the jet. 
The problem is formulated in three dimensions in order to investigate for- 
ward and backward scattering phenomena (7) , as well as the influence of iet 
flew speed. 


II FORMULATION .AND SOLUTION OF THE WDF.L 


The mechanisms by which shielding occurs are reflection and refraction 
of sound at the boundary between the jet and the surrounding air and by 
diffraction around the jet. 

The noise source is modelled by a stationary , discrete frequency point 
soiree located at ( r 0 » 9 0 > z 0 ) • The shielding jet is a cylinder of radius 
a, and is infinite in extent a-x.njt the z-a-.is. The temperature and flew ve- 
locity are uniform across the cylinder cross-section. Tne model is illus- 
trated in Figure 1. The expression for acoustic velocity potential is 
written for two regions; region I is outside the jet, region II is within 
the jet. 

In region I (outside the flew) 
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In region II ( inside the flow) 
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Where : 


(r , 6 q , z q ,) - location of point source 
u - source frequency 
q - source strength 
o - fluid density 
c - sound speed 

M - mach number - (iet flow speed/c ^) 
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Note: The subscript o refers to conditions outside the flow (ambient), and 1 
refers to conditions within the heated jet. 

The boundary conditions at the interface between the ambient air and the 
jet are: 

1) Pressure cont' : lity 

(P) 0 * (P) x at r * a 

-or - 

[■^Vo ■ 31 

2) Continuity of the vortex sheet (2). This condition states that the 
displacement of the medium is continuous and symmetrical at the boundary; 
r = a. Denoting this displacement by r = (z, t) , then: 
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at r = a 


- or - 

(n t ) 0 m + Vn z^l’ at r = a 

Eliminating time from equations 1 and 2 by assuming: 

<t>(r, 0, z, t) * Hr, 0 . z)e" lut 


Then: 
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With boundary conditions 

I) (ib»o 0 ^)o = ^Vl* * p i v *zh at r = a 


4) 


5) 

6 ) 


j 


d 


II) with similar assumption that 


n(z, t) - n '(2)e’ iajt 


(-iwn 1 ) * (-iun* ♦ Vnp ^ at r 3 a 


Using the Fourier transform: 

<j> 
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with corresponding inverse: 
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the equations become: 

I) (outside flow) 
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II) (inside flow) 

= 0 

with boundary conditions 

(i«P *) 3 (iPj'Ku - Vk z ))^ at r 3 a 
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The solution to the two-dimensional wave equation with a source at coordi- 
nates r , 0 Q is given in chapter 7 of Morse f, Ingard (5) . Thus, the trans- 
formed acoustic potential of the wave incident on the How cylinder is: 




iq o e 


‘ ik - z o 
^ o 


in 


Swp 


Z e m cos m (e-6 ) 
m o 

m«o 


/ vww 

I VVWo' 


r <• r 

o 

131 

r ” r o 


where: 
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r * distance from the z-axis to the receiver 
r Q = distance from the z-axis to the source 

In addition to the wave emitted from the source, there is a wave 
scattered fiom the flow cylinder, of the form: 
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where : 

Hj^fh^r) - Hankel function of the first kind = + iY m> which is chosen to 

ensure that the wave is outgoing. The wave transmitted into tie flew cyl- 
inder is the solution of the homogeneous two dimensional wave equation, and 
is of the form: 
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The form, J m , is chosen since the solution must be finite at the origin of 
the coordinate system. 

Applying the boundary conditions at r ■ a, where a < r Q 

* 

1) Pressure continuity 
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2) Continuity of the vortex sheet 


O- )• 

vy in' r 


+ t> ) 

vy sc'r 


k c 

,° 0 fii, ) 

k^C| ^trans ; r 


at r * a 


Then: 
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The primes denote differentiation of the Bessel function with respect 
to the argument. 

The transformed acoustic potential in the far field, at r > r , is: 
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where: 
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applying the inverse trans form, and inserting the tine dependence, the 
acoustic potential is: 
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SUTMARY OF ANALYSIS TO DATE 


1) Numerical Integration 

A digital computer program has been written to integrate equation 19 
numerically. The integration technique employed is a point Gaussian 
Quadrature. As a test, the technique has been used to evaluate the integral: 
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where : 

M 3 o th order Hankel function 
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This expression is the integral component of the solution for a point 
source located at the origin, radiating into free-space; where the receiver 
is located at the (cylindrical) coord inates,(r,e, z ) . This expression is 
chosen because: 

1) The integral is a constant for any combination of r and z, such that 

2 2 

■\. r + z * constant 


2) The magnitude of the integral decreases by one-half for each doubling 
of r 2 + z 2 . 

The Harikel function becomes unbounded as the argument approaches zero, 
when k z ■ k Q . Because of this discontinuity, the numerical integration 
scheme does not meet convergence criteria. Ongoing efforts in this area 
are concentrated in modifying the numerical integration technique to accom- 
modate the integration near a discontinuity. 

2^ Simplification to TVo - Dimcasions 

A further simplified, yet still significant model is a line source 
impinging on the cylinder of flow. In this model, the problem is reduced 
to two dimensions and the integral form of the solution is eliminated. The 
method of the solution is the same as for the three dimensional model. The 
velocity potential expressions are: 

Region I (exits ide flow) 
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Region II (inside flow) 
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where: 
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with boundary conditions: 

1) Pressure continuity 

■°o ( *t ) o " - p l ( Vl at r " a 

2) Normal velocity continuity 

(♦ r ) 0 ■ (* r ) ! at r « a 
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The wave incident upon the cylinder is 
-io)t 
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The wave scattered from the cylinder is: 

♦sc * e ' lUt i'\ c ° s( ”"' ,l1 ni < k o r ’ 
m»o 

The wave transmitted into the evlinder is: 
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The solution for the total (incident ♦ scattered) velocity potential in the 
far field (r > r ) is: 
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scattered 


As a check, when >> o Q (flow cylinder approaches a .solid cylinder) and 
r 0 > > (incident wave approaches a plane wave), the expression for ’■he 
scattered wave reduces to the solution for a plane wave incident upon a 
cylinder as given by Morse and Ingard (Sect. 8.1 of referenc (H ). 

Figure 2 is a plot of the directivity index based on equation 24. The 
directivity is 

PI(r s) * If) 1 oo , n 
o i in 



where: 

Tp,3 - coordinates of the receiver with respect to the roise source 
P-j.^ - incident plus scattered sound pressure from equation 24. 

D 

REF - sound pressure at the location without the How cv Under, the incident 
part of equation 24. 

The figure is drawn for the following conditions: 


11 



a) Outside flow - ambient conditions 

b) Inside flow cylinder - T * 2100° R, and Cj from ideal gas relation- 
ships . 

c) Jet diameter - 2a = 2.5 feet 

d) r Q - Jet separation - 5.0 feet 

e) r - radial distance from source - 40.0 feet to 1000 feet 

P 

f) Source frequency - 337. H . 

4 

The source frequency is evaluated based on a Strouhal Wimber of 0.25 for 
a Mach 1.5 heated jet, where the dominant source frequency is: 

M-c.- St 

f = ± 

d 


where: 

M - mach number = 


c^ = sound speed at 2100° R = 2250 ft. sec. 

St = Strouhal number 
d = jet diameter 

The directivity index is plotted for distances from the noise source of 
40, 125, 250, 1000. The contours show the barrier effect within the shadow 
tone, defined from ray acoustics as occuring between 165.5 and 194.5 degrees. 
The scattered wave is concentrated into side lobes distributed around the 
source. The lobe at 60 degrees corresponds to the critical angle from 
Snell's law, at which pure reflection is expected to occur. 

One interesting result is that the acoustic intensity is lower directly 
behind the noise source (e=o) than it is within the shadow zone. This 
indicates that little backscattering occurs at large angles of incidence, 
and that the sound wave passes through the jet. 


Future work with the two - dimensional model will he directed toward 
comparing the results obtained by the model to experimental and other 
analytical results reported in the literature (7,8). 
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SIAMARY 


A two - dimensional model of jet shielding has been developed and is 
currently being tested against other such models and experimental results. 
The purpose of development of the two - dimensional model has been to 
verify the process of the derivation and to provide insight into the devel- 
opment of the three - dimensional model. 

An equation for the acoustic velocity potential in three-space has 
been developed. Hie expression is in integral form. Current efforts have 
been centered on development of a numerical integration scheme to solve the 
expression. These efforts have been unsuccessful to date, but work continue 
in this area. An approximate closed form solution, incorporating a form of 
saddle point method is also being applied to the problem. This method has 
been successfully applied to radiation into free space from a point source. 
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